function V1 = Vec1(A, Ar, B, Br, varargin)
% Vec1 函数用于求解线性方程组 AV + IVAr + SUM_NVNr + BBr' = 0 中的未知变量向量 V
varargin1 = varargin(3);
n = size(A, 1); % 获取矩阵 A 的行数
r = size(Ar, 1); % 获取矩阵 Ar 的行数
result = zeros(n * r); % 初始化一个大小为 n*r 的零矩阵

result = result + kron(varargin{1}, varargin{2});
% 根据输入的 varargin 求和构建 result 矩阵
for i = 1:2:length(varargin1)
    result = result + kron(varargin1{1}{i,1}, varargin1{1}{i+1,1});
end

% 加上额外的两项，构建完整的线性方程组
result = result + kron(A, eye(size(Ar))) + kron(eye(size(A)), Ar);
B1 = -B * Br';
b = B1(:); % 计算方程组右侧的向量 b

v = result \ b; % 求解线性方程组，得到未知变量向量 v

V1 = reshape(v, n, r); % 将求解得到的 v 向量 reshape 成 n × r 的矩阵 V1
end
